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Abstract 

This paper presents a covariance analysis of the 
performance and sensitivity of the j*V**u C 
determination Extended Kalman Filter (EKF) used by the 
On Board Computer (OBC) of the Extreme Ultra Violet 
Explorer (EUVE) spacecraft. The linearized dynamics and 
measurement equations of the error states are derived 
which constitute the "truth model" describing the real 
behavior of the systems involved. The "design model 
used by the OBC EKF is then obtained by reducing the 
order of the truth model. The covariance matrix of the 
EKF which uses the reduced order model is not the 
correct covariance of the EKF estimation error. A true 
covariance analysis" has to be carried out In order to 
evaluate the correct accuracy of the OBC generated 
estimates. The results of such analysis are presented 
which indicate both the performance and the sensitivity 
of the OBC EKF. 


1.0 INTRODUCTION 

The Extreme Ultraviolet Explorer (EUVE) Is scheduled to 
be launched by a Delta launch vehicle in August 1990 
into a 550 km orbit with a 28.5 degree Inclination. The 
EUVE experiment will observe stellar objects emitting 
electromagnetic radiation with wavelengths of 100 to 
1000 angstroms. The spacecraft design is called an 

Explorer Platform (EP). The EP is designed tc ' be 
flexible enough to be used by many different 
experiments. The EP consists of three main modules: 
exDeriment, Platform Equipment Deck (PED), and the 
multimission modular spacecraft (MMS) . The MMS contains 
the attitude control system, power system, and the 
connand and data handling system. The EUVE mission is 
divided into two phases: all sky survey and 

spectroscopy. In the all sky survey, the spacecraft 
will be rotating at 3 revolutions per orbit (3 RPO 
about the roll axis) while instruments perpendicular to 
the roll axis scan the sky. Six months later, EUVE 
will be three axis stabilized at selected spectroscopic 
targets. 


2.0 ALGORITHM 


2.1 INTRODUCTION 

The attitude of the Explorer Platform (EP) is 
determined by gyros which measure the angular rate 
vector of the EP, by two fixed-head star trackers 
(FHST), and by one fine sun sensor (FSS). The gyros 
yield three components of the angular rate vector of 
the EP rotation with respect to inertial space. The 
components measured by the gyros are the projections of 
the vector on the body axes, which are the axes of the 
attitude control system (ACS). 
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If the exact orientation of the ACS with respect to 
inertial coordinates is known at some point, if the 
gyro outputs are perfect, and if no computation errors 
are introduced when solving the attitude propagation 
equations, then from that time on the EP attitude is 
known exactly. However, since the initial knowledge of 
the attitude is never perfect, since the gyro outputs 
include measurement and misalignment errors, and since 
the computation is not perfect either, the attitude of 
the EP is not perfectly known. Moreover, the attitude 
errors tend to diverge and consequently corrections of 
the computed attitude have to be performed. This is the 
reason for employing two FHSTs and one FSS in attitude 
determi nation. 

Star tracker and Sun sensor measurements, when used 
correctly, check the attitude error growth. The 
information supplied by the FHSTs and by the FSS is 
blended with the attitude computed based on the gyro 
outputs and on the initial EP orientation. This 
blending is done by a Kalman filter (KF). The EP on 
board computer (OBC) software uses quaternions for 
attitude determination. The relationship between vector 
measurements, which are the outputs of the FHSTs and of 
the FSS, and the quaternion of rotation is non-linear. 
Therefore an Extended Kalman filter (EKF), rather than 
a KF, has to be employed. 

For simplicity of implementation the EKF used by the 
OBC is actually a reduced order suboptimal filter which 
does not contain all the error sources in the gyros, 
In the FHSTs and in the FSS [1,2]. It is, therefore, 
necessary to investigate the predicted performance of 
the on board reduced order EKF. To evaluate the 
performance of the EKF, finer error models have to be 
used in describing the performance of the true hardware 
[3,4]. Such models are referred to in the literature as 
"truth models" [5]. In this paper we introduce a "truth 
model" which takes in account factors neglected in the 
OBC EKF model. These factors are gyro, FHST and FSS 
misalignments, gyro scale factor errors and the effect 
of the Sun not being captured in a narrow field of view 
about the boresight of the FSS. 

A convenient analysis tool is the "true covariance 
simulation [6]. This paper presents such analysis of 
the performance of the on board attitude determination 
EKF of the EUVE satellite. The "truth model" of the 
attitude determination problem is developed next. The 
"design model" is then listed in Section 2.3. Next the 
"true covariance* simulation algorithm is presented in 
Section 3. The analysis which was carried out and its 
results are presented in Section 4 and finally, the 
conclusions drawn from this analysis are presented in 
Section 5. 


2.2 IH£ TRUTH MODEL 
Error Propagation Model 

Open-loop attitude determination 

Consider Fig. 1 which describes a generic attitude 
control spacecraft (S/C). The input quaternion, q 1t 
represents a command attitude and q^. is the quaternion 
which represents the actual attitude of the S/C. As 
shown in Fig. 1, gyros which are mounted on the S/C 
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measure Its angular velocity. The readings of these 
gyros are used by the Attitude Determination algorithm 
to compute the quaternion q c which represents the 
computed attitude. This configuration is called "open- 
loop" since here q c is not fed back 


quaternions as depicted in Fig. 2. As indicated in this 
diagram 


q t = < 

H -to-a 

(la) 

q c - c 

H -to-c 

(l.b) 

dq " q a-to-c 

(l.c) 




J ATTITUDE CONTROLLED | q t 

| SPACECRAFT | 

i i 



I I 

| ATTITUDE | 
| DETERMINATION | 


<*c 

Fig. 1: Generic open-loop attitude determination 
configuration 


The subscript "t" corresponds to the transformation to 
the true attitude of the vehicle whereas the subscript 
M c" denotes the transformation to the computed 
attitude. (Note that H c" corresponds to the 
subscripted notation ’measured’ in refs. 2 and 3). When 
defining a quaternion of rotation and especially when 
dealing with quaternion products, a special care has to 
be given to the question of what coordinate frame the 
quaternion is referred to. If each of the three 
quaternions defined is referred to the coordinate 
system from which it transforms vectors, then the 
following relation between them holds 

<i-to-c * ^j-to-a ^a-to-c (2) 

where the product on the right hand side of (2) Is the 
quaternion product (defined in the Appendix) and the 
superscript^ denote the frame to which each quaternion 
is referred . From (l.c) the rightmost quaternion in 
(2) is dq. Note that dq is the only quaternion 
referred to the body frame whereas the other two are 
referred to the inertial frame. Keeping this in mind we 
omit all superscripts and use the notations of (1) to 
write (2) as 


into the attitude controlled S/C. Consequently in this 
configuration the gyro error, dw, of these gyros is not 
affecting the S/C attitude. 

Let us denote by "1* the inertial coordinate system 
which is the reference coordinate system and by N a" 
the ACS coordinate system which which we assume to be 
identical to the body system. The attitude 
determination problem is that of finding the quaternion 
which corresponds to the transformation matrix from "1" 
to "a” (or vice-versa). Since the gyros introduce 
measurement errors (dw) the computed attitude is 
erroneous. Therefore the computed transformation matrix 
which is supposed to transform vectors from the M i M 
to the ”a" frame, actually transforms the vectors 
from the "i" frame to another erroneous coordinate 
system which we denote by "c". Thus we distinguish 
between three coordinate systems; namely, the "1 ", the 
"a" and the "c" systems. We assume that the error in 
computing the transformation matrix is small, 
consequently "a" and "c" are almost identical. In 
other words, a very small transformation takes us from 





Differentiation of (3) yields 


q c = q t dq + q t dq (4) 


It is well known [7] that 


<U ’ l <lt w (5) 

where W is a quaternion of angular velocities defined 
as follows 


W - iw x + jWy + kw z (6) 


The components w x ,w y ,w z are the components of the true 
angular velocity vector at which the ACS coordinate 
frame rotates with respect to the inertial frame, 
coordi nati zed in the ACS frame. These components are 
measured by gyros which supply measured (and hence 
erroneous) data. In the lack of knowledge of the true 
rates, the gyro outputs are used in computing the 
quaternion, therefore the solution of (5) yields q 
rather than q t ; that is, q c is the solution of C 


Fig. 2: Schematic description of the quaternion 
relations in the attitude determination 
probl em 


the "a N to the "c" system. To the transformations 
between the coordinates there correspond suitable 


<3c 


2^c w m 


( 7 ) 


Note that q]_ t0 _ c - qj to . c qJ 
all quaternions are referred t v W1 
the order of the product is reverse! 


q ioHi 


; that is, when 
inertial frame 
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where W m is of the form of (6) only that the quaternion 
components are the measured rather than the exact 
angular rates. 

When (5) and (7) are substituted into (4) we obtain 

^c w m * ^ q t w dq + q t dq (8) 

Define a quaternion of angular-rate error as follows 

dW - w m - W 

then 

W - w m - dW (9) 

Substituting (9) and (3) into (8) yields 

^t dqW m * J q t( W n - dW > dq + q t dq 

which can be written as 

q t [dq + [(W m - dW)dq - ^dqWJ * 0 


d4 * -w^xdi + ^dedw + ^dwxd| (11* a) 
do = -dw .da (11* b) 


When dq expresses a small rotation, its vector 
part, dfc, is small, therefore the last term on the 
right hand side of (11. a) is of second order and hence 
is negligible. The right hand side of (11. b) is 
negligible too and indeed, since the absolute value of 
any quaternion of rotation is equal to 1, the scalar 
part of dq satisfies the equation de - [1 - |d£| ] ' 
and since the vector part is small, de, stays close to 
1, hence its time derivative is nearly zero. Note that 
as da is nearly 1, the second term on the right hand 
side of (11. a) is not negligible. Consequently (11) 
yields 

di ~ -w^xd* + ^dw (12. a) 


and as explained above 


de~l (12- b) 


The equation of interest is (12. a). 


Since q t is invertible, it is possible to pre-multlply 
both sides of the last equation by the inverse of q^. 
This yields the result that the expression In tne 
brackets is equal to zero and consequently 

| dq - ^dqV m - ^W m dq + |dWdq (10) 


Let us express the quaternions appearing in (10) in 
a more explicit form by their vector and scalar parts. 
Accordingly 


dq 


1 d| 

|.... 

I 

j_de_ 



where d| is the vector part of the quaternion. When the 
quaternion product is carried out (see the Appendix for 
the rules of quaternion product), (10) reads as follows 


The transformation matrix T a which corresponds to dq 
can be expressed in terms of small Euler angles. Define 
the angles as follows 

f is the roll angle error defined about the body x-axis 
u is the pitch angle error defined about the body y-axls 
y/ is the yaw angle error defined about the body z-axis. 

Note that for small rotations the order of rotation is 
irrelevant and we may refer all angles to the Initial 
coordinate system which prevailed before the small 
rotations took place. Also note that these angles are 
referred to the body frame, "a", as is implied in (2) 
and (3). When the transformation matrix from the body 
to the computed frame is expressed as a function of the 
three Euler angles defined above and when the angles 
approach zero, the transformation matrix becomes 


i 1 

J 

r 

-yy j 

| 

T c - j -r 

i 

r 

\ 


-r 

ij 



r 

i| 

cfcxw^ + 

i 

d0 !6!mj 

i! 

y^xdi + 

de *mj 

2 


1 

2 j 



I 

-d6.w 



-w.di 

_ 1 


j j dwxdi + dedw 


L -dw.di J 

The last quaternion equation is equivalent to two 
equations, one for the vector part of q and one for 
its scalar part. Using the following rules of vector 
product, AxB - -BxA and A.B = B.A, the two equations 
can be written as 


On the other hand, in terms of the components of dq, 
the upper right elements of T a are [8] 


tj _2 ' 2dqidq 2 + 2dq 3 dq 4 

(14. a) 

t l,3 * 2dq l dq 3 ' 2dq 2 dq 4 

(14. b) 

*2,1 ’ 2dq 2 dq 3 + 2dq l dq 4 

(14. c) 


The first term on the right hand side of each one of 
the above equations is of second order and hence is 
negligible. On the other hand, dq* appearing in the 
second term is the scalar part of dq which we denoted 
by da. As noted earlier this component is nearly equal 
to one. For these reasons (14) can be approximated as 
fol lows 
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t 1.2 

- zd <l3 

(15-a) 

‘1.3 

_ -2dq 2 

(15. b) 

‘2,3 

- 2d qi 

(15. c) 


Comparing (15) to the corresponding elements In (13) 
yields 


-*<p- 


|S/C ATTITUDE 
_J COMMAND 


Torque 


GENERATOR 


J S/C 
[DYNAMICS j_ 


Jr 


w 



d( »i * l V 

dq 2 * 


! -c 




◄- 


ATTITUDE l| 
DETERMINATION T 
I 


dq 3 - 

The above three components of q constitute the 
elements of thus 


Fig. 3: Generic rate command control loop part of a 
spacecraft attitude control system 


Consequently, w, the actual angular rate of the S/C 
differs from the desired angular rate vector, , by 
the gyro drift rate vector, dw; that is ’ 


I V> I 
ll I 

di » jl > 1 (16) 

Lyj 


and (12. a) can 

be written 

as 




iyi 
1 ! 

I 0 

I 

W Z 

V 

1 y 1 

I 1 

1 dw x | 


1 i 

| v | - 

1 1 

! -* z 

0 

w x 

1 1 

1 v 1 

i 1 

+ 1 dW y | 

(17) 

l_f_l 

L w y 

" w x 

0 _ 

1 I 

urn 

1 1 
l_d«zJ 



Closed-loop attitude detenni nation 

In the case of a closed-loop attitude determination, 
the S/C is maintained at a desired (possibly time- 

varying) attitude by a closed control loop which uses 
the gyro outputs to keep track of the S/C angular rate. 
This is shown in Fig. 3 in which a part of the control 
loop of a generic attitude control system of a S/C is 
presented. The purpose of this control loop is to force 
the S/C to follow a prescribed angular velocity vector, 
Wi, and in particular to maintain a constant attitude 
when w^=0. (Normally the commanded rate is a 

function of the difference between a commanded 
quaternion and the computed quaternion, q c ). We note 
from Fig. 3 that 

w “ w i - w c (18) 

The control loop is designed to force w to vanish, then 
*c “ -i (19) 

and since 

w c - w + dw (20) 

therefore 

w i - w + dw 


w - w.j - dw (21) 

and in particular when - 0 ; that is, when the S/C 
is required to maintain a constant attitude 

w - -dw (22) 


that is; the S/C drifts at the drift rate of the gyros 
but in a direction opposite to the gyro drift. We 
conclude from this discussion that * differs from the 
commanded rate by dw. Consequently the attitude of the 
S/C differs from the commanded attitude by the attitude 
error angles y, and y. In this case therefore, the 
attitude errors develop according to the following 
equation rather than according to (17) 


~° w z “y j j y j M"x| 

-w 2 O -x | | (23) 

_ w y -w x 0j|jf_| |_dw z _| 

Indeed when the commanded angular rate is zero (23) 
yields 


f= -dw x , -dw y and -dw z 

The right most term in (23) is not a white noise 
vector, therefore this dynamic model, while correct, is 
not suitable for use in a KF algorithm. To solve this 
difficulty the standard procedure of "signal shaping" 
is applied. This is done by considering the non-white 
vector as an output of a linear system whose input is 
white [5]. This is accomplished as follows. 

The elements dw x ,dw y ,dw z are the errors in measuring 
w. In other words, they are the errors in the x,y and 
z gyros respectively. It is assumed that an accurate 
enough model of the gyro errors is a one where there 
are five contributions to dw, which we denote by Mi » 
ds, ^ and Hi; that is, 1 

da - Mi + M 2 + + m (24. a) 

where Uj is a vector of constant drift rates of the 
gyros, M 2 a vector of random walk components of the 
gyros, Is the vector of gyro scale factor errors and 


1 >1 - 
I ! 

i_r_i 
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5 # 


is the error due to gyro misalignments, nj the 
ite noise component of the gyros. Since Uj is 
constant 

1 1 - 0 (24. b) 

and since is random walk then we can write 

1 2 - Il2 (24. c) 

where a? is white noise. Note that despite the 
notationf n? does not have the units of angular 

velocity. Let us denote the vector of scale factor 

errors by k, then 


& * * [k x , ^y* 


(24. d) 


where * denotes the transpose and k x , ky and k z are the 
scale factor errors of the x, / and z gyros 
respectively. The expression for di is given by 


! w x 


0 I 


di - (0 w v 0 | X 

I ^ I 

jo 0 W 2 j 

Note that since K is a constant 

k • o 


(24. e) 


(24. f) 


The gyro errors due to misalignments are generated by 
the projection on the gyro Input axis of the angular 
velocity components which are nominally perpendicular 
to that axis. The misalignment angles are the angles by 
which the gyro sensitive axis is off from its nominal 
orthogonal direction towards the other two coordinate 
axes. Consequently we have 


j w jix 

- J^y 
j w *z 


|w v m 7 0 0 0 0 

I y 

|0 0 w z w x 0 0 


10 0 0 0 


I l?xy 

11/ 

II 


xz 


*x "yj |*yz 

j^yx 

j*zx 

!^zy 


(24. g) 


of the gyro outputs in order to obtain more accurate 
gyro measurements, therefore the same result is 
achieved when we estimate their sum. When we combine uj 
and into one state denoted by n, we may use (24. b) 
and ( Z4. c) to write the dynamic model of y 


U - C2 


(24.1) 


In order to augment the models presented in (23) and 
(24) we define the following matrices 


! * ! 

1 0 

w z 

-w y 

-1 

0 

0 


i i 
1 » i 


0 

w x 

0 

-1 

0 


i i 

1 r ! 

I 

,* i w y 

" w x 

0 

0 

0 

-1 


i i 

A - | 







! u * ! 

1 0 
i 

0 

0 

0 

0 

0 


! u > ! 

1 o 

0 

0 

0 

0 

0 


L u z-I 

L o 

0 

0 

0 

0 

0 



(25) 


-n * 


n lx 



0 0 

Wy- 

■W Z 0 

0 

0 

o 1 
1 

"ly 


i o 

-Wy 0 

0 

0 -W Z 

-*X 

0 

o 1 
1 

"lz 


A** - | 

0 -W Z 

0 

0 0 

0 

- W X 

-Wyj 

n 2x 


1 o 

0 0 

0 

0 0 

0 

0 

o 1 
1 

"2y 


I o 
1 

0 0 

0 

0 0 

0 

0 

o 1 

1 

- n 2z- 


l_0 

0 0 

0 

0 0 

0 

0 

0J 







... 


(26) 

t' 

m ^xy» 

**xz* ^yz* 

^yx’ t 

zx» 

*zy] 



(27) 


We may also want to consider the misalignment angles of 
the two FHSTs and of the FSS, therefore let us denote 
the vector of. the three misalignment angles of the 
first FHST by 1 < > that of the second FHST by c t and the 
vector of the FSS misalignment angles by where 


h- 

- C'^X- 



(28. a) 

h' 

■ [ 2 *x* 

h r 


(28. b) 

s i’ 

* 



(28. c) 


where 1-a,y.z J-x,y,z Is the non -orthogonality 
angle between the i-th gyro and the j-th axis. Now 
since is constant we can write 

(24. h) 


*1J 


1-x,y,z 

j-x,y,z 


The next step in the derivation of the dynamics matrix 
is the augmentation of the system error model given In 
(23) with the gyro error model given in (24) [5] . bucn 
an augmented model has, in our case, 27 states. 
Fortunately, we can eliminate 3 states by combining the 
constant drift rate components and the random walk 
components into one error. This will th * 
possibility of distinguishing between them, but this is 
of no great consequence since even If we can estimate 
them separately, we subtract them both from the reading 


Since all of these angles are constants we may write 

1* . 0 h • 0 s i - 0 (29) 


With the above information and notations we can now 
write the augmented dynamics equation of the "truth 
model \ The augmented dynamics equation of the "truth 
model" Is given in (30). The validity of (30) can be 
verified by examining (23) - (29). 
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l“ 1 1 


1 X 1 
1 1 

A* 

A** 

0 

0 

0 

* 1 

II X 1 

1 +11 
1 n II 


! 1 

i a i 
| | 

0 

0 

0 

0 

0 

1 1 1 
II X 1 

i 1 1 

12 11 
i 1 1 

jd 

1 I 

i * i 

0 

0 

0 

0 

0 

1 1 1 
II * 1 

1 1 1 
2 II 

|dt 

1 a 

i 1 * i 

1 l 

0 

0 

0 

0 

0 

111 1 + 

II 1 * 1 
1 1 1 

II 

1 2 II 


1 9 1 

i 2 * i 

0 

0 

0 

0 

0 

1 1 5 I 
II 2 * 1 

i i 

1 1 1 
1 2 II 

t f 1 


1 c 1 
I s * i 

Is 

0 

0 

0 

0 

0 

lie ' 
II s * 1 

i 1 1 

1 2 II 


1 1 

- 





J L J 

L J 1 


.... (30) 

Measurement Model 


Star tracer measurements 


* -Si-(^st)a + e y (32.b) 

We use an under-bar and a subscript to denote a column 
matrix whose elements are the components of the vector 
in question when resolved in the coordinate system 
denoted by the subscript. 

The observables (also known as effective 
measurements) which are processed by the EKF are the 
difference between measured and computed quantities. 
For star tracker measurements we feed the EKF with the 
difference between the measured components x and y 
given in (32) and the corresponding computed values 
which are obtained by transforming the star vector from 
inertial to body coordinates. The star vector 
coordi nat i zed in the inertial frame, which we denote by 
Sj, is precisely known from the almanac. We do not 
know, however, the exact value of T‘, the 
transformation matrix from Inertial to body system. All 
we know is the computed transformation matrix T*. The 
relationship between the two matrices is given by 


Define a star tracker coordinate system as shown in 
Fig. 4. The z axis points along the boresight of the 
star tracker. Consequently the x and y axes ajfe in 
the image-plane of the star tracker. Denote by 5 the 
vector In the direction of the star and whose length is 
the length of the light path from the image-plane to 
the optics. It is assumed that the light which Is 
emitted from the star towards which the star tracker is 
pointing, hits the image-plane close to the boresight 
such that it can be assumed that the distance between 
the optics and the image plane Is nearly equal to that 
of the light path from the optics to the image-plane, 
i.e. h~|S|. The signals measured by the tracker are the 
projections of on * 
the x and the y 
axes which, as 

mentioned, are in 
the Image-plane. The 
OBC converts the 

two outputs of the 
star tracker to 
tangents of A and B. 

Obviously, the tan- 
gents of A and of B 
are, respectively, 
the projections of 
on the x and y 
axes of the tracker, 
where 5 is a unit 

vector in the 
direction of the 

star. That is, if we 
denote these axes by 

— st and y s fc 
respectively, then 


-£.x 


st + e x 
... (31. a) 

-S-*st + e y 
... (31. b) 






1 c 1 a 


(33) 

where 

Tq is the 

error 

matrix 

given in (13). 

If we 

define 

the matrix 0 

as follows 




1 o 

r 

-»-i 



0 ■ 

1 . f 

0 

Y 1 

I 

(34) 



L* 

-J* 

o_l 


then 


T? 

i + e 


(35) 


We also do not know the exact direction of jj s t and 
since the star tracker Is misaligned. We only know 

Xst.ass. and *st,ass. which are the vectors and y st 
In the assumed coordinate system of the FHST (that 
Is, in the non misaligned FHST). Consequently, the 
computed values are calculated In correspondence with 
(32) as follows 


x c * U s t,ass)a 

(36. a) 

*c * *(Tc$i) . (ist,ass^a 

(36. b) 

Using {33} and (35) these two equations become 

x c - -[(I+0)T*Si].(x st , ass ) a 

(37. a) 

Vc * ■ [ (I+8)TjSi ) . (y s t t ass^a 

(37. b) 

We note that 


T i*l - S a 

(38) 

therefore (37) can be written as 


x c “ '^a + 0 £a3*(*st,ass)a 

(39. a) 

+ 6 £al • (*st,ass)a 

(39. b) 

When we now difference (32) and (39) the 
equations are obtained 

following 


where x m and y m are the tangents of A and B, " ,H 
denotes the dot product of vectors and e x and 
are measurement noise signals which are assumed to be 
zero-mean white processes. (Actually, the OBC converts 
x ffl and y m into components of the unit vector £. For 
small A and B these components are basically equal to 
the respective tangents). Let us now express the vector 
quantities of (31) In the EP body coordinate system 

x n * -Sa-l*it>a + e x < 32 - a > 


Z 1 " - *c 

z 2 * - *c 

We note that 


"^a'(*st)a + e x + + • (X s t t ass^a 

... (40. a) 

■-^a’^st^a + e y + tS* + 0£ a ] * (^st,ass^a 

... (40. b) 
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*st - M J'‘” -st.ass 


*st 


J.ass 

m|’ ass 


*st,ass 


(41. a) 
(41. b) 


where M? ,ass Is the transformation matrix from the 
assumed FHST coordinate system to the actual one. In 


analogy to (35) it can also be shown that 
Mj.ass . , 


where 


- [*x] 


I 0 
I 

! -h 


{&] 


(42) 


(43) 


is the transformation matrix from the assumrd tracker 
frame to the body frame. This matrix is known 
precisely. Let us write 

hJ.«s - [ ujj, e 2 , m 3 ] (50) 

where mi, ffl 2 *"d m 3 are the three columns of Mj- a “. 

It is clear that (X st ,ass>T- the unit ve . ct0 . r n 
tracker assumed x axis expressed in tracker 

coordinates, is given by 


(*st,ass)T 


I 1 I 
I 0 | 
I 0 I 


0 

L ty 'K 

The angles t\ i*x»y>z are the misalignment angles of 
the actual FHST coordinates with respect to the assumed 
FHST coordinates. Note that because of the closeness of 
the two, the angles are the same in either coordinates. 
(See development leading to (56)). When we substitute 


and similarly 


I 0 | 

(*st,ass)T ■ jJJ 


(51. a) 


(51. b) 


Therefore when (50) and (51) are substituted into (49), 
the following is obtained 


(42) into (41) we obtain 


(*st,a$s)a * ®1 

(52. a) 

x s t “ ( I ’ l^ x ) ) ^st.ass 

(44. a) 

^st,ass^a * E 2 

(52. b) 

y s t ■ ( 1 " (i)x] ) ^st.ass 

(44. b) 

When (521 is substituted Into (46), and when 

S a is 

hence 


replaced in (46) by $ c which is given in (4B), tne 


following is obtained 


*st " 2 st,ass ‘ & x ^st, ass 

(45. a) 

Zi - Sc-lMa x ill + eSc-El + e x 

(53. a) 

Yst ’ list, ass ' x y st,ass 

(45. b) 

z 2 - Sc-lW), x m 2 ] + 8Sc ®2 + e y 

(53. b) 

When (45) are substituted into (40) the following is 
obtained 

Noting that 


z x - -Vtist.ass ‘ i X x st,ass>a + e x + IS, + 


X 

id 

OH 

« 

CD 


♦ eSjMXst, 

ass^a 

and using the vector Identity 


z 2 - -S a -(y s t,ass ' 4 x ^st.ass^a + e y + tS a + 


(A xB).C» (B x C).A 


+ * Cx$ t , 

ass^a 

equations (53) can be written as 


which, after some multiplications and subtractions, 

zj - (n^ x S c )^ a + x -c) *®a + e x 

(54. a) 

yields 

(46. a) 

z 2 » (m 2 x S c ).^ a + (rn 2 * £ c )*®a + e y 

(54. b) 

Z\ * £*♦(£ x ^st,ass^a + ®^a -^st,ass^a + e x 

Expressed by its components t given as follows 


z 2 “ 5 a *(i x *st,as$)a + e ^a- ^st ,ass^a * e y 

(46. b) 

r*fi 

i * i 



For (46) to be useful, we need to evaluate S a , 
<X st ass>a and ^t,as$)a- As mentioned earlier, 

Is Mown from the a Tma’nac, therefore 

Sa - Tl S, « 7 > 


I 

I dy 
I y 
U>z- 


(55) 


We do not know Tj but we do know T c 


whicji, for 

small attitude errors, is quite close to T a . 


The 


dllid I I OIL i twuc * - “ ,1 , , a r 

replacement of the true value by its estimate is one of 
the features of an EKF, so, we too, follow this 
practice, compute 


however, the angles which constitute the components of 
are defined in the FHST coordinate system. Ther|fore, 
to keep using the same angles we writ e s d a - M T r , , A 
But M| ,ass * I - [£x], therefore Mj- 1 i m & lC x J6 
Consequently, 


T{ Si 


(48) 


and use £ c rather than S a in (46). Next we handle 
the computation of (Xst.assU ancl (^st.assla - ' s 
clear that 

(49. a) 
(49. b) 


(*st,ass)a “ M a ,aSS ( x st,ass)T,ass 
(^st.ass^a ’ ^a’ aSS ( Jt st,assJT,ass 


i a - * 

When (56) is substituted into (54) we obtain 
zj - (mi x S C M + (mi x ScJ.fij + « x 
z 2 . ( E2 x S c ).4 + (m 2 x Sc). 0 a + e y 

Me can write (57) as follows 


(56) 

(57. a) 
(57. b) 


where T denotes the star tracker coordinate «st|m 
defined at the beginning of this section and M a ’ 
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f(®i x $ c r ii y> i 

111 + 

Ullfe x S C )’_||U | 

I I 
L'f_l 


I (Ej x 5^) 1 1 dx I I e x 

+ I II x l + I 

J_(m 2 X £c>’J I * X I I e y 

\J Z \ l_0 _ 


... (58. a) 

where * denotes the transpose. This equation is the 
measurement model of a generic FHST. Each one of the 
two FHST has such a measurement equation. Writing (58) 
in terms of the state vector X and thus forming the 
measurement matrix, H, Is straight forward. The matrix 
which corresponds to this measurement equation for the 
first FHSS is 


i I (mi x $ c )* 

s - 1 

l_<n»2 x 2 C )’ 


I (mi x SJ’ | -| 

°2X 1 2 I I °2x6 I 

I (n>2 x Sc)’ I _| 


... (58. b) 


and the matrix for the second FHST is 


, I (E] x S c )’ I 

Z «F “ I I 


M Z x Sc)’ 


(mi x Sc)’ I "I 

, C w ! °2x3 ] 

(m 2 x | j 

... (58. c) 


The geometry of the FSS measurement is similar to 
that of the FHST presented in Fig. 4. Here, however, we 
cannot assume that the angles A and 8 are small; that 
Is, the Sun vector Is not nearly coincidental with the 
boresight line. Therefore all the developments that 
were based on this assumption are not valid in the 
development of the FSS error model. Consequently a 
different approach has to be taken. It Is evident that 
£ shown in Fig. 4, can be expressed In the FSS 
coordinates, s, as follows 

$5 - [-tanA, -tanB, Ijd (59. a) 

d - [(tanA) 2 + (tanB) 2 + l]' 1 / 2 (59. b) 

Let £ s ass denote a column matrix whose elements are 
the components of £ in the assumed (non-misal igned) FSS 
coordinates. The relationship between this vector and 
5 S is given by 

$s • G !’ aSS Ss,ass (60) 

where G|* ass is the transformation matrix from the 
assumed to the fine Sun sensor coordinates. In analogy 
to (42) we can write 

G|’ ass - I - [ s *x] (61) 


% -%\ 


Substitution of (61) into (60) yields 

* £ s>ass - [ S l)x] S s>ass ( 63 ) 

From (63) we immediately realize that if instead of 
£$,as$ we use Ss> we introduce an error due to the FSS 
misalignment. This error is - [ s dxl S 

1 * J -s , ass * 

The FSS outputs are really (tanA) m and (tanB) m where 
(see fig. 4) tanA ■ S x /h tanB - S v /h. The subscript 
m denotes the measured tanA and lanB. Define the 
following column matrix 

S;,m- (-( tanA ) m > "(tanB) m , l]d m (64. a) 

d m - [(tanA)* + (tanB)* + l]' 1 /* (64. b) 

L *l.> ‘ (tanA) m and v m - (tanB)_, then (64) can be 
written as 


£s,m“ £' u m» ~ v m> ^ d m 
d m * (“m + *m + I]' 1 /* 


Furthermore, using Taylor series expansion we can write 
S s,m ■ S s + d 5 (66) 



"“11 

W 12 

‘ir v 
1 



- “21 

<M 

CM 

31 

l 1 

l'-v 

- «£ 


-“31 

“23- 

f 


W Is 

evaluated 

as 

follows 

[9] 

“11 - 

d m ' d m u 

S 

“12 * 

VVm 

“21 * 

“l2 


W 22 * 

d m " d m v m 

“31 ■ 

-4 


“32 ■ 

' d m v m 


and e x e y ® r ® the additive measurement errors 
Involved in melisurlng tanA and tanB respectively. From 
(66) and (67) we obtain 

*s,m “ $s + w * (69) 

Substitution of (63) into (69) yields 

-s,m " £$,ass ’ £$, a $s + Wf (70) 

Next we compute the estimate of We denote the 
computed value by £ s c . The computation is carried out 
as follows 

•^s,c * G s,ass T c^I (71) 

In (71) we actually have to use the matrix T ^ • 
however, since this matrix is unknown to us we use a T J 
instead. £j is taken from the ephemeris. From (33) - (35^ 

T c - {> - tfix])T| 

Thus (71) can be written as 

S s,C - G|, ass {I - [6x]}T's r 


which can be written as 


*s,c - “s.assSa * G f.ass^xji a 

where £ a Is a column matrix whose elements are 
the components of $ when the latter Is resolved In the 
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body, a, coordinate system. Note that since Gf jas s is 
orthogonal 

- ^ iass [6x]S a - - (Gf >ass fi) x G| iass S a - 

- - (6S.*ssfl) x -s,as 


Note that the 3rd row of the square matrices in (77) is 
omitted. This completes the development of the "truth 
model". To sum it up, the dynamics model 1$ given by 
(30), the FHST measurement model is given by (58) and 
it fits either one of the two FHSTs, and finally, the 
measurement model of the FSS is given by (77). 


Using the last equation, we can write the former as 

S s ,c " ^a.ass ■ ( G s,ass£l x -s,ass ( 72) 

Now define the first component of the effective 
measurement (which is to be processed by the EKF) as 
follows 

(73) 


Z 1 


^s,m " -s,c^x 


where { } Y denotes the x component of the expression in 
the brackets. When (70) and (72) are substituted into 
(73) , we obtain 


Z 1 " ass ’ -s, ass + 

- 1 


which yields 
Z 1 


a, ass + ^ G s,ass 


£) x ^s,ass^x 


<- [Ss.ass^lassS + Ss,ass x + W §>x < 74 ) 


Following the rationale that led to (56), we can 
substitute S s ass in (74) by S s _ with practically no 
loss of accuracy. In addition, use the notation 

G s,ass ' l3l» 3 2- 33] (75) 

therefore (74) can be written as 

*1 - U-Ss,m x][3i. 32. a 3 ]0 + S s ,m x H + we) x (76. a) 

To compute z ? we apply the foregoing development but 
now we use the y rather than the x component. This 
will yield a result similar to (76. a); namely, 

*2 - u-£ s , m xJIai. 32- 3 3 ]9 + Ss.m x + Wfi) y (76 ‘ b) 


r 1 

! z 2 


i 

I 

l3ixS s>m . 32xS S(rn . S 3 xS s J 
-1 2 




I. 


I 

Jz 


\y\ 


|v! + 
1 | 


Wj 


W 

L 1 + w 2 

, IN 

r _ ii 

INI 

INI 

H 

1 0 II 
1 1 


1 1 1 


... (77. a) 


2.3. HE DESIGN MODEL 
Error Propagation Model 

The “design model" is the simplified, reduced order 
model which is assumed to be the model of the system 
for the OBC filter design purposes. The following 
assumptions are made in the design of the EUVE Update 
Filter. The gyro scale factor errors and misalignments 
are negligible (or fully compensated for). The FHSTs 
and the FSS are perfectly aligned. With these 
assumptions the error propagation equation of the 
"truth model" reduces to 



0 

0 

If 

1 

>"! 

I w ^ x I 

1 

0 

1 1 
II 

7t- | 

I 

j w ly j 

0 

-1 

1 1 
II 
II 


"lz j 

0 

0 

II 
1 1 

u x | 

| w 2x | 

0 

0 

I I 

II 
1 1 

u y j 

w 2y j 

0 

0 


- U z-I 

L w 2z-l 


... (78) 


This model can be expressed as 

l* - aY + n* (79) 

where X*, A* and n* are defined in (25). 


Measurement Model 

Star tracker measurements 

With the assumption mentioned before and 
corresponding to the "design model" of (78), The FHST 
measurement matrix of either star tracker reduces to 


I r-f] . T'" ^ >T ! !>!.!"< ! 

j |_ Z 2 _I l _( rn 2 X S c ):i|V t I e y | 

i i_y_! i_o j 

This yields the following measurement matrix 


(80. a) 


where the subscript 2 denotes the first two rows of a 
matrix. The measurement matrix, H s which corresponds to 
this measurement equation is given by 


I (ffij x S c )’ | 

L ®2 x S c )’l 


°2x3 


(80. b) 


r i 


1 1 
1 1 

It is easily seen that this is the measurement 
for either FHST. 

matrix 

H s - l9l x ^s,m’ 32 x ^s,m’ 93 x £$ ( nij 

i 

°3x24 

j t£ s> m x l j 

I 1 o 

Sun sensor measurements 


L ' 


1 - 1 2 
... (77. b) 

Corresponding to the state vector of the 
model" the FSS measurement equation is reduced 

"design 

to 
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... (81. a) 

The corresponding measurement matrix is 

, f I "I 

H s- laixSs.m- a 2 xS Stm , g 3 x$ s J 0 3x3 | (8i.b) 

i- I J 2 

3.0 1RU£ COVARIANCE SIMULATION ALGORITHM 

To present the so-called "true covariance" 
simulation algorithm whose development is introduced in 
[6], we have to define, D, the transformation matrix 
from the state vector of the "truth model" to that of 
the "design model". It is easy to see that in our case 

D • f'exel °6xl8 ) (82) 

Using 0 we define the following matrices 

dA - DA - A*D (83. a) 

where A is tjje 24x24 matrix defined in (30). The 
matrixes A, A and dA are then used to define A c as 
follows 


i A | 0 | 
| — | | 
_dA| A | 


(83. b) 


Next, we discretize the matrices A c and Q where the 
latter is the spectral density matrix of the white 
noise vector driving the dynamics part of the "design 
model" given in (78). The discretization algorithm is 
given in [5, pp. 296 - 299]. The discretization is 
denoted as follows 


A C — > Cl 

Q* — > Cl 

With Q k _j on hand, we compute Qj'.j as follows 


Q Li 


i i i 
i— i 
I D ( 


l"i I d t I 

L I J 


(84) 


In our case, all the preceding matrices are constant 
and need to be computed only once. 

Between measurement updates we propagate the 
matrices C k and P k as follows 

c k(-> - »i-i c k-i *tr + Ci ( 85 - a ) 

p kH - *k-l p k-l( + ) *k-l ’ + Qk-l (85. b) 

where C k is the second moment matrix of the augmented 
state vector whose entries are, from top to bottom, the 
state vector of the "truth model" given in (30), and a 
vector which is the difference between the state 
eslmate generated by the OBC EKF and the correct value 
of this state. Note that the second vector, which has 6 
components, is the correct estimation error vector. 
Therefore the last 6 elements on the main diagonal of 
C k are the mean squre errors of the filter estimation 
error and their evaluation is the goal of the "true 
covariance" analysis. In contrast to these 6 elements, 


the 6 elements on the main diagonal of P k are the 
apparent variances of the estimation error states. That 
Is, If the "truth model" were identical to the "design 
model", these elements would have been the variances of 
the estimation error. 


When a measurement is acquired, the following 
computations are carried out 

K k s p k(-) H k* t H k p k(-K’ * ^r 1 (86. a) 

p k (+) - [I-kJhJ] Pj(-) [I-kJhJ]’ + K k R k K k ’ (86. b) 

dH k = H k - H k D {86. c) 


C k (+) - B k C k (-) B k ’ + < k R k (86. e) 


4.0 CASE SIUPY OBJECTIVES AND RESULTS 

There were three primary objectives in the case 
studies. First, the performance of the EP filter was 
examined in the ideal situation when its model was 
equivalent to the truth model. Secondly, the expected 
onorbit behavior of the filter was examined. Lastly, a 
sensitivity analysis was performed. The cases studied 
were as follows: 

- Case 1: No Errors 

- Case 2: Expected Errors 

- Case 3: Sensitivity Analysis 

- 3A: Gyro white noise about each axis 

* 3B: Gyro random walk about each axis 

- 3C: Gyro Misalignments about each axis 

- 30: Gyro Scale factor errors about each axis 

- 3E: FHST noise 

- 3F: FHST # 2 misalignments about each axis 

• 3G: FSS noise 

- 3H: FSS misalignments about each axis 

Each simulation was ten minutes. The attitude and 
gyro drift estimation errors were determined by the 
truth model and update filter. The results in sections 
4.1, 4.2, and 4.3 represent the truth model determined 
estimation errors. The following is a listing of 
nominal simulation Input values: 

I niti al State Variances 

Initial Attitude Error: 1800 arcsec/axis 

Initial Drift Rate Bias: 0.5 arcsec/sec/axi s 

Dynamic Noise Inputs 


IRU White Noise drift (roll): (0.68936 arcsec/sec 

f.ji.i i _ . , : • 


-r - 


1/2 2 


(pitch and yaw axes): (4.246E-2 arcsec/sec 1 ^)* 

IRU Random Walk Drift: (4.4413E-5 arcsec/sec 5 / V/axis 


Measurement Noise Input 


FHST Measurment Noise Variances: (14 arcsec) 2 

FSS Measurement Noise Variances: (24.4131 arcsec) 2 


4.1 Case J.: No errors 

The following case demonstates the performance of 
the filter in the ideal case when the truth model was 
identical to the design model. The final attitude and 
gyro drift errors were as follows: 
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Attitude Estimation Errors 
(arcseconds) 

Roll Pitch Yaw 
4.8784 2.8321 2.5080 


Gyro Drift Estimation Errors 

(arcseconds/second) 


Roll 


Pitch 
n QimF-9 


Yaw 

0 01*P_ 


4.2 Case Z: Expected ElCflli 

This case demonstrates the expected performance of 
the filter on orbit. The attitude is defined relative 
to one of the FHSTs (#1 in our simulations). Thus, the 
obtained attitude accuracy is on the order of the 
accuracy of the FHSTs which are the primary attitude 
sensor. The following are additional expected on orbit 
input errors: 


3C x 

3C y 

3C z 

3D x 

30 y 

3D z 

3E 

3F x 

3F y 

3F Z 

3G 

3H x 

3H y 

3H z 


16.5005 

16.5005 

16.5005 

16.5005 

16.5005 

16.5005 

19.1600 

27.2160 

22.3209 

17.3819 

16.9075 

16.7568 

16.8615 

16.5005 


19.1984 

19.1984 

19.1984 

19.1984 

19.1984 

19.1984 

25.7578 

29.6328 

25.5485 

19.3448 

18.9557 

24.9939 

19.8543 

19.1984 


13.0119 

13.0119 

13.0119 

13.0119 

13.0119 

13.0119 

18.0298 

15.9332 

13.0534 

21.7191 

13.0284 

13.1199 

16.3459 

13.0119 


FHST #1 Misalignment: 0 arcseconds/axi s 
FHST #2 Misalignment: 24 arcseconds/axi s 
FSS Misalignment: 36 arcseconds/axi s 

Gyro Scale Factor Error: 1000 ppm/axis 

Gyro Misalignment: 8 arcseconds/axis 

The results are the following: 


Attitude Estimation Errors 
(arcseconds ) 

Roll Pitch Yaw 

16.5005 19.1984 13.0119 

Gyro Drift Estimation Errors 
(arcseconds/second) 

Roll Pitch Yaw 

3 . 2513E- 2 5.1160E-2 6.9354E-2 


Comparing these values to the ideal case, one can see 
errors Induced by only considering a subset of the true 
state vector in state estimation. The attitude 
estimation errors are off by several orders of 
magnitude and the gyro drift estimation errors are off 
almost an order of magnitude. 


4.3 Case J: Sensitivity Analy sis. 

In the following simulations, the sensitivity of 
onboard filter to additional attitude sensor noises, 
misalignments, and scale factor errors was tested. 
These errors were applied separately to each sensor 
axis, and the resulting attitude and gyro drift 
estimation errors were observed. Sensitivity to the 
various error sources were determined in the following 
manner about each spacecraft axis (where applicable): 

Dynamic noise (white & random walk): 3x nominal/axis 

FHST 02 Misalignments: 2x nomi nal/axi s 

FSS Misalignments: 2x nominal/axis 

Measurement noise (FHSTs & FSS): 2x nominal 



00 300 00 ' 

TIME (SECONDS) 


■iq. 5: The standard deviation of the pitch estimation 

error vs. time. (The solid line is that of the 
true error, and the dashed line is of the 
error predicted by the update filter 
covari ance matrix. ) 

The time history of the standard deviations of the 
attitude and gyro drift estimation errors was plotted 
for Case 2 (expected onorbit errors . The results 
about each axis were found to be similar. A typica 
nlot (Ditch attitude error) is presented in Fig. b. 


The following tables and figures demonstrate the 
filter performance due to the increased errors. 




Attitude 

Estimation Errors ( ercsecondil 

Case 

Axis 

Roll 

Pitch 

Yaw 

3A 

x 

18.6886 

19.1984 

13.0119 

3A 

y 

16.5006 

19.2332 

13.0179 

3A 

Z 

16.5005 

19.2065 

13 .0435 

38 

X 

16.5007 

19.1984 

13.0119 

3B 

v 

16.5005 

19.1994 

13.0121 

3B 

J 

z 

16.5005 

19.1987 

13.0125 


Gvro Drift Est1"9tl£n EIT2TS 
(arcseconds/ second) 


Case 

Axis 

Roll 

Pitch 

Yaw 

3A 

3A 

3A 

X 

y 

z 

8.7164E-2 

3.2513E-2 

3.2513E-2 

5.1160E-2 

5.1451E-2 

5.1202E-2 

6.9354E-2 

6.9364E-2 

6.9566E-2 

3B 

3B 

3B 

X 

y 

z 

3.2564E-2 

3.2513E-2 

3.2513E-2 

5.U60E-2 
5. 1197E-2 
5.1161E-2 

6.9354E-2 

6.9355E-2 

6.9378E-2 
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3C 

X 

3 . 2513E-2 

3C 

y 

3. 2513E-2 

3C 

z 

3.2513E-2 

3D 

X 

3.8577E-2 

3D 

y 

3 . 2513E-2 

3D 

z 

3 . 2513E-2 

3E 


3 . 5639E-2 

3F 

X 

3.2513E-2 

3F 

y 

3 . 2513E-2 

3F 

z 

3.2513E-2 


5.1160E-2 

6.9354E-2 

6.8993E-2 

6.9354E-2 

5. 1160E-2 

8.3385E-2 

5. 1160E-2 

6.9354E-2 

5. 1160E-2 

6.9354E-2 

5. 1160E-2 

6.9354E-2 

6.5980E-2 

8.9894E-2 

5.9675E-2 

1 0 . 2473E-2 

5. 1277E-2 

8.9347E-2 

7.7423E-2 

6. 9809E-2 



Fig. 6: The standard deviation of the pitch gyro drift 

estimation error vs. time. {The solid line is 
that of the true error, and the dashed is of 
the error predicted by the update filter 
covariance matrix.) 


3G 3 * 2543 E - 2 

3H x 3.2513E-2 

3H y 3. 2513E-2 

3H z 3 . 2513E-2 

A typical plot for gyro 
Fig. 6. 


5. 1259E-2 6.8545E-2 

5. 1467E-2 8.7576E-2 

6.0912E-2 7 . 1389E-2 

5.1160E-2 6.9354E-2 

drift error is presented in 


Of the gyro noises, the white noise component had the 
most effect on the filter performance. As expected, the 
effect was confined primarily to the axis being 
corrupted. When the X-axis gyro white noise was 
increased by 3x, the roll gyro drift estimation error 
jumped 5.0E-2 arcseconds/second, and the roll 
estimation errors jumped approximately 2 arcseconds. 
The effect of an Increase In pitch and yaw gyro white 

?? 1se !\ ad * ver * nom1nal *ffect on estimation errors. 
Since EUVE has a roll rate of 3 RPO, the roll gyro has 
to use less accurate gyro data as compared to the pitch 
and yaw axes which are approximately Inertial. This 
inaccuracy in roll gyro data Is modeled by an increase 
in the white noise component In the roll gyro data 
Thus, an increase In the white noise about the roll 
axis affects the attitude much more significantly than 
an increase about the pitch and yaw axes. The gyro 
misalignments about the pitch and yaw axes corrupted 


the r respective drift estimates significantly due to 
p ™ Ject1 1 ° 1 ns P. 1ck1n 9 up the relatively high rol? 
rate. The roll axis misalignments have no effect due 
to zero yaw and pitch rates. The gyro scale factors 

up in .^ he i; 011 9yn> drift estimation error 
due the above mentioned high relavlve roll rate and 0 

the li?" r a , te .V. The FHST n "^alignments were 
the largest contrlbuter to attitude and gyro drift 

estimation errors as expected with the FHST X-axis 
misalignment causing roll and pitch errors of 27.2 and 
29.6 arcseconds. The FHST Z-axIs caused a yaw error of 

a i t Th€ , FS ?_ mfsal ignments affect on the 

attitude and gyro drift estimation errors were 

thp" /Jr 3nt bUt n0t * S s 1 9 n 1f1cant as the FHSTs due to 
the larger sensor noise variance. The FSS X-axis 

misalignment translated Into a pitch error of 24.9 
vau Se ptt?. *!, hne 3 r ' ax 1 s misalignment caused roll and 

respectively^ 'fhTruu 16 8 and 16 ' 3 ^seconds 
respectively. The FHST measurement noise increases 

as the Fi^c? at . ti ^ de estimation errors almost as much 

errors of lJ .""’l 1 i 9 n " ,e " ts 10 w * th Pitch, and yaw 

rors of 19,1, 25.7, and 18.0 arcseconds respectively 
An increase in FSS measurement noise had a relatively 

lull eff H Ct 0n es A? a V 0n errors - The resulting roW 

pitch, and yaw attitude errors were 16.9, 18 9 and 

estL a Hnp eC ° ndS re ? p * ct1vel y- The pitch attitude 
estimation error went down slightly as compared to the 

nominal simulation, and the roll and yaw attitude 

inCreaSed S,l9ht1y ’ The reas °" the attitude 
estimation errors were affected so little as compared 
to the increased FHST noise simulation was due to the 

heavier Than'Vh'^ l thS "' 0re accurate FHST measurements 
heavier than the less accurate FSS measurements. 


5.0 CONCLUSIONS 

r ° F VJthe errors, the FHST misalignments proved to 
o.M 6 7i be most significant attitude and gyro drift 

CZTh The : o11 and p,tch estimation 
errors increased by approximately 10 arcseconds from 

the nominal estimation errors when the X-axis 
misalignment of FHST #2 was doubled to 48 arcseconds 

ya 9 W h e 8 stl'mat1on 1Sa1i9n, " ent K 0f FHST n lncreised 
me yaw estimation errors by approximately 8 

e r rr S or C s 0n were IffeK a "\ ya "A° 'dMftenVmVtion 

errors were affected most by 1 and Y - axis 
misalignments of FHST #2. An increase in the whit! 

?a™~t ab0U li the 9y r° X - axis " as responsible for hi 
r ? 11 ..^ 0 dr1ft estimation error. The gyrl 
drift estimation errors only affects the system when 

r^f e iX t UPdat K e Per l 0dS are lar 9 e - EUVE should have 

the FHStI If nU n m nt er tl updates from 

rn?,iH H !i T !i f ® t ’ the attltude estimation accuracy 
could degrade significantly. Overall, the results 
showed the EUVE update filter to be quite robust evel 

svstem h S Thf<r S . 9 H 1f J Cant errors were Put into the 
in Up' fin study demonstrated the six states modeled 
in the filter are the most significant states needed 
for onboard attitude estimation 


Appendix 


Express the two quaternions qj and q 2 as follows 


9l 


I I 
I I] I 
!■■■ I 
I Pi I 


and 


^2 


Q2 

r "i 

I r 2 | 

I... I 
I P 2 1 

l:i 


where £j, r, are the 
quaternions and pj, p 2 


vector parts of the respective 
are their scalar parts. Then 
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I 

I Il x E2 + V\tz + P2El 


qiP2 


P1P2 


clr 2 


The upper part of the column yields three components 
which are the components of the imaginary part of the 
quaternion product and the lower part yields the scalar 
part of the product. 
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